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ABSTRACT 


Th«  linear  theory  of  elaaticity  ia  uaed  to  inveatigata  axially  aymmetric  wave 
propagation  in  an  infinitely  long  two-layered  cylinder.  Each  material  ia 
taken  to  be  homogeneou*  and  iaotropic.  A  perfect  bond  ia  aaaumed  at  the 
interface^  while  the  inner  and  outer  boundaries  of  the  composite  cylinder 
are  treated  as  traction-free.  The  diepereioo  determinant  relating  phase 
velocity  and  wave  number  for  a  harmonic  train  of  waves  satfslying  these 
boundary  ..auditions  is  presented.  The  character  of  the  dispersion  equation 
is  investigated  aisalytically  and  numerically.  Stress  and  displacement 
distributions  are  also  prei  ented  for  the  numerical  example.  Comparisons 
arc  made  with  sui  approximate  solution  of  the  same  problem  obtained  by 
means  of  a  thin  cbeli  theory  incorporating  thicknese -shear  daformation  of 
each  layer. 
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1.  INTRODUCTION 


In  recnnt  ynnra  consldarnbl*  nttnntlon  has  baan  focuaad  on  multi- 
layarad  sheila.  Often  a  sandwich-type  construction  is  used  to  lighten  tha 
weight  of  a  shell  structure;  in  other  instances,  a  protective  layer  is  bonded 
to  a  shell  as,  for  example,  in  the  case  of  a  re-entry  vehicle  heat  shield  or 
a  rocket  noasle  liner.  Two-layered  Uddi  shell  configurations  are  also  used 
in  solid  propellant  rocket  motors. 

This  study  employs  the  linear  theory  of  elasticity  and  treats  the  propa¬ 
gation  of  a  train  of  waves  in  an  infinitely  long,  two-layered  cyliitder;  each 
layer  ia  homogeneous  and  iaotropic.  Computationa  Ijased  on  this  solution 
are  uaed  In  an  assesament  of  tha  accuracy  of  an  approximate  two-laycrad 
shell  theory  that  wae  preeented  recently.^  Thia  latter  the<<ry  included  the 
effects  of  shear  deformation  and  rotatory  inertia.  In  Ref.  1,  a  partial 
comparieon  is  made  between  the  shell  theory  and  the  solutions  obteined  in 
ths  present  work  ueing  the  linear  theory  of  elaaticity.  Dispsraion  curvas 
were  compared  in  detail,  and  a  few  displacement  distribution  comparisons 
were  made.  The  present  work  givec  detailed  displacement  comps risons 
over  e  larger  frequency  and  wave  nomber  regime  ar.d  also  presents  strsss 
distribution  comparisons. 

The  propegalion  of  waves  in  cylindrically  bounded  media  has  been 
extensively  Inveetigeted.  Althoa|^  most  of  the  work  has  been  limited  to 
cylinders  of  a  single  meterial,  it  ia  of  interest  lo  recall  a  few  of  the  more 
pertinent  references.  Pochhammer^  and  Chree^  first  formulated  the  prob¬ 
lem  for  solid  cylindrical  hers.  Ghosh^  formulaiert  the  problem  (or  hollow 

cylindrical  bars  but  preaented  no  calculations.  L^ter,  Gaxis^’  and 

e  q 

Gresntpon  *  made  extensive  numerical  celcutations  for  the  vibretinns  of  e 


For  e  more  extensive  bibliography  dealing  with  the  thin  shell  literature  of 
layered  sheila  aae  Ref.  I. 


hollow  cylinder  and  compared  them  with  veverai  approximate  ahell  theories 
■  iirK  *■  tKoB*  nf  Marrmanm  an<(  lulirakv^^'  arwt  Naahdi  an<l  Cjwtomr .  ^ 

Other  inveatigatora  include  Bird^^  and  Bird,  Hart,  and  McClure, 

Toe  eibrafiona  of  a  multilayered  cylinder  using  the  equations  ot  the  lineal 

theory  of  elaaticity  have  not  been  extensively  investigated.  Baltnikonis, 

17 

Gottenberg,  and  Schreiner  treated  simple  thickness -shear  vibrations  of  a 
two-layered  cylinder,  and  McNivsn,  Sackman,  and  Shah^^  treated  propagation 
of  axially  ayanmctric  waves  in  solid  bars  with  an  outer  finite  layer.  Other  than 
these,  the  authors  know  of  no  other  references  treating  multilayered  cylinders 
by  the  linear  theory  of  elasticity. 

19 

Layered  half-planes  have  been  extensively  treated  by  geophysicists, 
but  their  work  is  not  particularly  of  intereet  hare  since  geophysical  earth 
modale  always  have  one  infinite  layer.  Of  the  plans  two  medium  problems 
with  finite  layers,  the  symmetrical  sandwich  two-dimensional  bsam  is  treated 
by  Salto  and  Sato.  and  the  asymmetrical  two-layerad  counterpart  is  treated 
hy  Jones.  The  lest  refsrencs  is  especially  apropos  sines  It  is  shown  hsrs 
that  the  wave  propagation  solution  of  the  two-leyered  cylinder  problem  degen¬ 
erates  into  the  solution  ot  the  plane  two  layered  medium  problem  when  the 
wavelength  becomee  sufficiently  email  compared  to  the  thickness. 

The  present  anelyeic  is  formulated  in  terms  of  displacsment  potentials. 

A  solution  in  ths  form  of  an  Infinits  train  of  axially  symmetric  waves  is  assumed. 
To  satisfy  boundary  conditions  the  {diase  velocity  (or  frequency)  must  depend  on 
wave  number  in  such  a  way  that  sn  eighth-order  determinant  vanishes.  Due  to 
its  complexity  little  analytical  progress  can  be  made  with  this  dispersion 
detsrmlaaat  excapt  in  special  caaas.  However,  for  infinitely  long  waves  the 
determinant  reduces  to  a  product  of  two  fourth- order  determinants  whose 
frequency  roots  correspond  to  vibratione  with  eithsr  pursly  axial  or  pursly 
radial  motion.  Alternatively,  for  very  short  waves  the  determinant  reduces  to 
a  form  given  previously  by  Jonee^^  for  a  plane  two- layered  medium.  Hers  the 
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phaae  velocity  root*  correspond  to  Rayleigh  waves  on  the  free  surfaces  and 
a  noaaihle  .Stnneley  wave  at  the  interface.  For  waves  of  intermediate  length, 
roots  of  the  dispersion  determinant  are  found  numerically  with  a  digital  com¬ 
puter  prugrarn  <a«ai.iiu«u  hersin.  Dlsplaccmnr.t  ar^  stress  distributions  cor¬ 
responding  to  these  roots  are  also  found  for  a  specific  numerical  example,  and 
these  are  used  to  estimate  the  accuracy  of  the  previously  mentioned  shell  theory. 

It  is  concluded  that  a  Timoshenko-type  shell  theory  gives  good  agreement 
with  the  present  exact  solution  in  a  region  of  applicability  encompaasing  low 
enough  frequencies  snd  'arge  enough  wavelengths.  It  appears  that  to  extend 
this  region  of  applicability  one  must  use  s  shell  theory  incorporating  thickness- 
stretch  niotio...  On  the  basis  of  the  displacements  obtained  from  the  exact 
theory,  a  linear  distribution  of  radial  motion  does  not  appear  to  be  an 
unreasonable  firat  approximation  for  a  thickness -at retch  theory. 


$ - 

A  Stonsiey  weve  mey  or  may  not  exist  at  an  interface  depending  on  the 
elastic  properties  and  densities  of  the  two  media  (Kef.  19,  p.  113). 
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n.  ANALYSIS 


I 


Consider  a  doubly  infinite  hollow  cylinder  cmnpoeed  of  two  homo¬ 
geneous,  isotropic,  elastic  media  with  Lame  constaate  Xj,  P|,  X.2t 
uid  densities  and  where  subscript  I  refers  to  the  outer  layer 
and  subscript  2  refers  to  the  inner  layer.  The  cylinders  are  perfectly 
bonded  together  at  the  interface.  Cylindrical  coordinates  r,  8,  end  a 
are  employed.  The  interface  radius  is  denoted  by  a,  and  the  thicknesses 
are  h^  and  h^. 

Written  in  terms  of  the  potential  functions  ^  and  ♦  the  equations  of 
elasticity  for  motions  with  torsionleea  axial  symmetry  are 

a  P 

where  the  displacements  and  stresses  may  be  generated  from  tbn  potential 
functions  by 
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The  tertris  u  and  w  are  the  radial  and  axial  dlaplacementa,  reapactivaly, 
and  <T^^,  tr^^,  and  are  the  atreaa  coznpoaeiita.  Tha  expraaaion  for 
ia  omitted  since  it  is  not  eaaantial  to  the  problem.  All  other  atreaa 
components  are  zero  and  a  and  3  ara  tha  dilatational  and  aaitivnlurnimal 
wave  apeeria  defined  by 


P 


(6) 


Equationa  <i)-(5)  hold  in  cither  layer  provided  appropriate  values  of  a, 

X  and  |j.  are  used. 

If  a  Rayleigh  train  of  waves  ia  assumed,  then  all  quantit  ?s  are  con 
sidered  to  vary  as 


.  i(kz-wt) 

w  =  w(r)«  '  '  , 

asnd  the  solution  to  Eqa.  (1)  ia  (the 
compactness) 


'rr  " 
factor 


(kS'Uit) 


etc. 


(T) 


has  been  suppressad  for 


f  =  AI^j(kcr)  +  BKjjfktr) 
ClQ(k6r)  4  DKg{k6r) 


<8) 


Here  and  are  modified  Beosel  functions  of  the  first  and  second  kind  and 
zeroth  order.  The  solutions  Eq.  (8)  hold  in  either  medium  provided  the 
parametere  appropriate  to  each  medium  are  uoed  in  the  equations.  Thus  there 
will  be  two  sets  of  c's  und  &'s: 
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(9) 

(cont.) 
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wh«re 
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Z 
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(10) 


Also  there  will  be  two  sets  of  constsnts  Ap  A2  And  Bj ,  B2>  to  be 

determined  by  application  of  boundary  conditions  on  o^^.u^and  w. 

FJxpressing  the  displacements  and  stresses  in  terms  of  the  solutions,  Eq.  (6), 
one  obtains 


1  +  62)i^(k,r)  -  2.2 


l,(k«r)1  r  2  2 

+  B  (1  +  6Vo(k*r)  +  2.^-^ 


2  2  r  IlOtOr)-}  r  K,(k6r)-l| 

+  2k262  C^Io(k6r)  -  j  .  +  -Wr-J 


^11 


(11) 


rz 

P 


2ik^.[AIj(k«r)  «  BK,(k«r)j  +  ik26[c(l  +  .  d(1  ♦  62)Kj(k6r)] 

(12) 

Uj.  =  k^A.IjOt.r)  -  B«Kj(k.r)  +  C6I,(k6r)  -  D6Kj(k6r)j  (13) 


w 


ikfAI^tk.r)  ♦  BKo(\.r)  4  C62l^(k6r)  4  D62KQ(k6r)]  (14) 


The  above  expressions  are  valid  in  each  region  provided  appropriate  values 
of  A. ,  B. ,  C, ,  D. ,  . . ,  6 . ,  u  .  (i  =  1 . 2)  are  used. 


Th«  boundary  conditions  for  tbs  £r«#  vibration  problam  ara 
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r  =  »  -  h. 


<r<^)  =  .<2) 
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TIm  aupBracript*  indicate  quantitiaa  in  medium  t  or  2.  From  the  boundary 
coodltioe,  Eq.  (15).  one  obtains  with  the  aid  of  Eqs.  (1 1)>(i4)  eight  linear 
hrnnogeneous  equations  for  the  eight  constants  Aj.  .  .  02-  Since  the  squatioiu 
are  long  and  their  formulstions  ere  strsighiforwsrd,  the  squstions  are 
presented  in  the  Appendix  rsthsr  than  in  the  text. 

To  assure  nontrivial  values  of  A|,  Bj,  .  ...  C2>  D2>  the  determinant 
of  their  coefficients  must  be  set  equal  to  ssro.  This  constitutes  the  dispersion 
equation.  The  determinant  is  as  follows; 


mi'* 


Tha  determinant,  Eq.  (16).  it  to  compticatod  that  Uttla  can  be  dotto  to 

—  ta  n  %  ff  lat  a  I  I  a  •»  Adaanaa%  aw  4  a  waa  a  m4  A4aA  iftA^aaoMHa  d 

nant  art  obtainable  for  the  limitinf  caact  of  very  long  and  very  abort  wavea. 

For  infinite  wavelengtht,  i.e.  ,  k  0.  the  tolutlona.  Sq.  (B).  bacotne 
independent  of  r,  and  the  determinant  dagenoratea  to  the  product  of  two 
fourth>order  determinants.  A  mo^e  straightforward  darivatiMit  of  theae 
fourth-order  determinants  involves  repeating  the  calculationa  witii  the  s 
dependence  excluded  from  the  outaet.  Then  it  ia  clear  that  one  of  tha  fourth- 
order  determinanta  corrasponda  to  purely  radial  motions  while  the  other 
corroeponde  to  purely  exial  motions. 

Equating  the  radisi  motion  determinant  to  aero  one  obtains  the  frequency 
equation  for  eimple  thickneaa  stretch  vibrations: 


Eqaatiaf  Ui«  ajdftl  motioii  to  Mro  OM  olittiM  Hm  tv^qmmuY 


aquation  for  axial 

abaar  vibratiaM 
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'‘i[^  '*  -  ‘■i’] 

(18) 

C4{aakioo  (18)  h&s  b««n  obUiiMd  by  BAitmkoaia,  •!  al. 

For  v«ry  abort  wovalongtha  cotnporo4  to  tho  latorfoca  rodiaa  of  th« 
eyllndar,  it  would  b«  oxpoctod  that  propagatloo  of  aalaUy  aymmatrle  wavaa  io 
a  two-layarad  cylindar  would  diffar  vary  llttla  from  propagatiofi  of  atraight 
craatad  aravea  in  a  plana  taro^layarad  madiuan  at  laaat  for  a  thin  cylindar. 

Tba  truth  of  tbia  auppoaition  can  ba  ahown  analytically  in  Eq.  (16)  by  ra> 
placing  Ig.  Ip  Kq,  and  Kj  by  thair  aaymptotic  valuaa  for  largo  argumanta: 


(19) 


(*  »  1) 
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With  these  substitutions  and  a  slight  redefinition  of  coefficients  in  Eqs.  (A-1)- 

21 

(A-8),  the  dispersion  determinant  becomes  the  same  as  that  given  by  Jones 
for  the  plane  two-layered  problem.  As  would  be  expected,  the  waves  degen¬ 
erate  into  two  Rayleigh  surface  waves,  one  in  each  medium,  and  a  possible 
Stoneley  wave  at  the  interface. 

For  wavelengths  of  intermediate  size  it  is  difficult  to  extract  much 
information  about  the  character  of  the  possible  wave  propagation  solutions 
except  by  considering  specific  numerical  examples.  This  is  done  in  the 
following  section  where  the  results  of  the  numerical  example  are  also  used 
to  check  the  accuracy  of  the  shell  theory  given  in  Ref.  1. 
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m.  CCMFUTATIOKS 


In  ordnir  to  invottigata  th«  charnctnr  ot  tlM  wnv«  prr>yn§ntion  jolttUona 
fovnrncd  by  Cqa.  (A-l)>(A-8)  numnrical  •aalyaia  waa  acnployari.  A  prograon 
for  tba  IBM  7094  coanpotar  waa  davalopad  at  tha  Aaroapaoa  Corporntio«i  for 
datarmining  phaaa  valocitiaa  (and  from  tham,  fraquanclaa)  aatlafying  Cq.  (I6). 
To  facilitate  tha  uaa  of  tha  praaaat  aolotion  tor  aaaaaaiof  tha  range  of  validity 
of  approidmata  ahall  thaoriaa,  diaplacamant  and  atraaa  diatributiona  through 
the  cylinder  thichneaa  are  calculated  as  well. 

For  a  given  numerical  problem  valid  computer  raaulta  are  obtained  only 
for  a  limited  range  of  wavelengths.  Tha  largest  aumbar  available  lor  rootina 
calculation  is  10^*  exp  (t8)»  while  for  large  wave  numbers  tha  Ta  (modified 
Bessel  functions  of  the  first  kind)  are  ot  the  order  of  asp  (  0tH)(a/H)) .  Hare 
K  >  h|  4  h2  is  the  total  wall  thickness  of  tha  cylinder.  Thus,  for  valid  com¬ 
puter  rasulto.  kH  cannot  be  much  larger  than  M(H/a).  In  the  example  to  be 
considered  hers  a/H  is  30  so  that  kH  <  5  gives  a  (air  estimate  of  the  range 
amenable  to  computations  based  on  Cq.  (16).  For  larger  values  of  kH  the 
wavelength  is  short  compared  to  the  radius  of  the  cylinder,  the  asymptotic 
expressions  <cJ  Cq.  (19)  are  appropriate,  and  satisfactory  numerical  results 
are  obtained  using  a  computer  program  based  on  the  plane  medium  equations 
of  Rsf.  21. 

properties  of  the  cylinder  chosen  for  detailed  numerical  study  are  given 
in  Table  I.  This  example  was  used  to  chock  the  validity  of  the  apprcndmats 
shall  thsory  of  Del.  1  whors  dispersion  curves  wsrv  compered  for  low 
frequencies  and  large  wavelengthe.  Additional  dispersion  curves  are  prassntsd 
hers  for  the  first  nlns  modss  of  sxislly  symmetric  wsvs  propagation.  Dis¬ 
tributions  of  displacamsnts  and  strsssss  through  As  thicknsss  of  ths  cylinder 
are  also  presented,  and  for  the  firet  four  modas  theee  are  compared  with  the 
dietributions  p^edirted  by  the  epproxtmete  ehcU  theory.  ^ 
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Tabl*  1.  Nuro*ric«l  Prcf«rtia*  at  tli*  C^ind«r  Studied 

t  TntiAr  Lav#:^ 

KUterial  Conatante  »  4.00  X  loSb/ia. lO^tb/la.^ 

Vj  *  1/4  =  1/3 

*»1  *  (tSt ^  <»2  * 

Geometrical  Parameters 

^  =30  -jj-  --  0.  3  if  =  0. 7  <H  =  h,  ♦  hj) 


A.  DlSPBtRSlON  CURVES 

Oiapereion  curves  for  the  cylinder  of  Table  1  appear  in  Flga.  1  and  2.  In 

both  graphs,  the  abscissa  is  nondiinenaionsl  wave  number  A  s  kH.  In  Fig.  1  the 

ordinate  is  nondimenaional  phet^e  velocity  a  *  c/fl*  >  w/kB*,  and  in  Fig.  2  it  ia 

e  e 

nondimenaional  frequency  R  s  wH/O  .  The  reference  velocity  p  (P|  't 
The  approach  of  the  phase  velocity  to  various  limiting  values  may  be 
examined  in  Fig.  1.  For  instance,  fot  vory  long  waves,  tlie  phase  velocity  cf 
the  first  mode  approximates  that  of  "bar"  waves.  Also,  for  very  short  waves 
the  phase  velocity  of  the  first  mode  approaches  the  Rayleigh  wave  epeed  of  the 

s 

elower  medium.  For  the  second  and  higher  modes  the  phase  velocity  becomes 
very  large  as  the  wave  number  is  made  smaller.  Therefore,  the  long  wave¬ 
length  behavior  of  these  modes  it  better  exantined  in  the  frequency  wave  number 
plot  of  Pig.  2. 

In  Fig.  2  it  is  observed  that  curves  for  all  the  modes  except  the  first  have 
a  finite  frequency  intercept  for  sero  wave  number.  These  frequencies  and  the 
character  of  the  motion  (or  each  of  the  modes  are  euminariscd  in  Table  II. 

Other  features  of  ths  curvss  prssented  in  Fig,  2  Includs  the  relative  minimum 


-16- 


lO  CM  —  o 

•  *Aiio(n3A  3$vHd  immmmoH 


Fill.  2.  Frequency  ve  Wave  Number 


-1»- 


Table  11.  Character  of  the  Infinite  Wavelength  Axially  Symmetric 
V'ibrati''>'ia  of  the  Taro*  Layered  Cylinder  of  Table  I 


n 

Mode  Frequence 


Character  oi  Motion 


2 

3 

4 

5 

6 

7 

8 
9 


nodea.  axial  motion 

C.0fa6  sere  nodes,  radial  motion,  "ring"  vibration 

3.86  1  node,  axial  motion,  thickness  shsar 

S.80  .2  node u,  axial  motion,  thieknsss  shear 

7,07  I  node,  radial  motion,  thicknesa  stretch 

10.2  )  n-'xles,  axial  motion,  thicknesa  shear 

1  1.  I  2  nodes,  radi.al  metion,  thickness  stretch 

n.  I  4  nodes,  axinl  motion,  thickness  sheer 

15,9  5  nodes,  axial  motion .  thicknaas  shear 


of  mode  6  near  A  =  2.5  (this  correaponde  to  esro  group  valocity)  and  the  close 
approach  of  modes  7  and  6  near  A  =  3.  I.  The  inaet  of  Fig.  2  shows  this  close 
approach  on  a  magnified  scale. 

B.  DlSPUliCEMENTS 

In  Ref.  1  an  approximate  shell  theory  ie  developed  and  parti* *ly  compared 
with  the  present  exact  solutions.  Figures  3  throt  ^h  7  extend  the  comparison  of 
displacements  over  the  entire  range  of  interest.  Due  to  the  eigonvalue  nature 
of  the  problem  coneidcred  the  absolute  magnitudes  of  the  displacenients  are 
undetermined  since  muitipUeetion  by  a  ronetant  (normalisation)  factor  is 
permissible.  Therefore,  our  comparisons  are  of  the  shapes  of  the  displace¬ 
ment  distributions,  the  magnitude  having  been  adjusted  to  make  tha  prcsenl 
exact  theory  and  the  shell  theory  of  Rsf.  i  Kg  css  at  some  convenient  value  of  r. 
Solution  for  both  theories  of  a  forced  motion  problem  would  permit  comparisons 
of  magnitudes  as  wall.  However,  examination  of  the  shapes  alone  of  the  dis¬ 
placement  distributions  is  of  interest  since  formulations  of  higher  order  shall 
theories  generally  postulate  functional  forms  for  the  dependence  of  the  dis¬ 
placements  on  the  thickness  coordinate. 


RADIAL  OISPLACEMERTS 


AXIAL  OISHACEMCNTS 


-SHEaTwemr  ~~  elasticity  TMEcmy 

rif.  4.  Saeoad  Moda  DtefUMmaat  Dlatflbutlea* 


For  £ho  pr«««at  compAtisotu  th«  udiU  dUpUcwnWRtt  ha^ht  kvm 
oonnaliiod  to  make  «iioct  UMtory  sad  ahaU  tkaory  d^BjtlacanaqaM  afaal  to 
on«  at  th«  outor  tuiiacc.  Tliii  normalisation,  of  •roiirao,  fixoa  <lafiafto 
rauoa  notwaan  too  raaiai  mapiacomontB  of  too  two  liaooriaa  «v«z-f 
location  in  tho  ahoU  (aseapt  in  tha  caao  of  isfinita  sravalaagth  adMra  aaial 
and  radial  motion  ara  uncouplad).  Tha  radial  diaplacamaato  naad  ao'idiara 
ba  equal.  Rather  than  praaaat  than*  in  auch  a  (on»  wa  h*va  takan  tha  libaity 
of  introdttcittc  a  diffarant  noitnaliaation  of  tha  radial  diaplacamonta .  Tha 
radial  diaplacamanta  from  both  thaoriaa  ara  made  equal  to  one  at  tha  outer 
aurface.  Such  a  noimalisatioa  ia  nacasaary  for  uncouplad  radial  motion, 
and  nioraover  it  permits  tha  uae  of  a  uniform  acala  for  nearly  all  tha  radial 
displacement  plots.  The  lnforn;ation  suppressed  by  this  normallaatioa 
convention,  that  is,  the  ratio  of  the  maximum  outer  surface  radial  displace* 
mant  to  maximum  outer  surface  axial  dieplacament,  la  praaantad  In  Tnhla  HI. 

In  Ref.  1  the  firat  mode  axial  and  radial  diaplaeamaata  jaa  woU  na  tha 
diaparalon  curvaa)  ware  compared  for  a  limited  raaga  of  A  <  3.0.  Figure  3 
extanda  thia  compariaon  up  to  A  a  9.0.  Since  several  curvaa  for  0<  A  <3.0 
wera  praaantad  in  Ref.  1.  the  curvea  in  Fig.  3  atari  at  A  ■  5.0.  Further, 
for  A  <  3.0,  there  ia  little  variation  of  radial  diaplaeamant  thrott(^  tha 
thlcknaaa.  and  tha  axial  diaplaeamant  dlifara  Uttla  from  the  U'liaanr  dia> 
tribution  aaaumad  in  tha  ahall  theory.  Starting  at  A  «  5. 0  the  noouidlortnity 
of  radial  diaplaeamant  beglna  to  become  pronounced,  and  the  axial  diaplace* 
mant  becomes  increasingly  nonlinear. 

Previous  analysaa^'  show  that  aa  tha  wave  number  bacoEnas  incraaa* 
ingly  large  tha  lira!  mode  dacaya  from  a  flexural  mode  into  a  Rayleigh  srava 
In  tha  alowar  medium.  It  ia  apparent  la  the  curvaa  for  A  «  9.0  toat  thia 
tranaitloa  ia  nearly  complete:  tha  motioa  is  oeskcsatratad  almost  aatiraly  ia 
tha  outer  medium.  One  might  expect  certain  of  the  dAaplacamaat  plata  to 
axhibic  other  auifaca  or  interface  ssavaa.  However,  with  tha  paramatara  used 
for  these  two  naadia.  ealculatieas  ehaw  that  Stonalay  wavaa  d«  not  exiat  at  tha 
interface  <Raf.  19.  p.  113).  Theory  pradieta  that  for  aufflciaadly  high  wave 
numbers  a  Rayleigh  wave  will  form  ia  the  faatar  (inaar)  medium.  Shch  a  wave 
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Table  m.  Batioa  of  Outer  Surfaeo  Ifaalinta 

IMaplaoumaa*  to  Outer  aatfaca  iiaxiM 

ru^^^ - 

B  Badtal 

Iteo ' 

hfoda 

Wavo  Numbor,  A 

Iteaat  Tboory 

note.ilMBo 

Ibotl  TbAory 

S 

l.OM 

1.474 

1 

7 

1.440 

1.841 

9 

1.4M 

1,739 

0 

• 

m 

2 

1 

-0. 277 

-0.0145 

3 

>2. 421 

-0,048 

0 

0 

0 

3 

1 

•0. 1024 

-0.M47 

3 

-0.0424 

•41.0421 

0 

0 

0 

4 

1 

0.  422 

0, 0277 

3 

1.844 

0,0942 

0 

*  a 

S 

1 

-11.94 

-- 

3 

•21.40 

•  • 

22 

it  aiMUcattd  la  t  mort  comflieatad  way  caUad  "tarractaf ,  aad  »ta  oea* 
■idaratioa  would  raqulra  iavaatlgatloa  of  a  largar  raaga  of  A  thaa  kao  boon 
uadaitakoa  for  tli«  proaoni  work. 

Figaro  4  ahowr  tba  oaccod  nodo  diaplacanooia  for  A  ■  0.  1.  aad  3. 

For  A  *>  0,  tho  motloo  io  totally  radial,  aad  ao  plot  of  axial  dioplacomaata  to 
Bocaooary,  la  cootraot  to  tho  firot  taodo,  Uio  dlaplacoinoiat  <vurvaa  dagoaorato 
from  tko  aliotl  tboory  quit*  rapidly  oritli  a  pkaao  roraroal  la  tho  radial  dlaplaco 
laaat  alraady  apFarowt  atAal.O.  ByAsl.O  tho  aoaliaoarlty  of  th«  earrao  it 
proBoaacad.  Thay  aro  totally  dilforom  from  thooo  of  tho  ohoU  thoory.  TUo 
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rmdi»l  dii»]>lac«tn«nt  carv«»  for  both  A  «  1.0  and  3.0  ara  avidanca  that  ttaa* 
•idarabla  thickoasa-atratch  dafennatien  ia  praeaid. 

Kigurva  7  and  o  oii9w  in«  oiaiMacvmvinv  loc  iav  inirv  ana  tvvfM 

Thay  ara  pradotniaantly  tbickaaaa'ahaar  modaa.  At  A  «  0  dM  motlan  ia 
axial^  tharafora  oo  radial  diaplacamant  plota  apftaar.  Ia  both  fifuraa.  tha 
pradictioaa  from  shall  theory  cotnpara  qulta  wall  arith  tha  axact  solotloat  at 
laaat  for  tha  udal  diaplacaniants.  As  in  tha  sacoad  mod#  dtara  is  aetfieiaat 
thicknass- stretch  motion  present  to  cause  tha  appearance  of  nodes  In  most  of 
the  radial  displacement  plots.  Note  in  Pig.  5  that  tha  radisd  dis^acemaats 
for  A  s  3  ara  plotted  to  a  different  scale  than  tha  rest. 

Figure  7  shows  the  displacements  for  the  fifth  mode,  tha  first  mode  not 
predicted  by  ehell  theory.  At  A  =  0  the  motion  is  purely  radial  or  thlchnasa- 
stratch.  At  least  for  small  wave  number  tha  ahapas  of  tha  radial  dlsplaeamant 
curves  ara  aoeh  that  a  linear  distribution  might  offer  a  good  approximation. 
Thua.  a  higher  order  shall  theory  admittinf  a  linear  distrlbiitioa  of  radial 
displacements  might  give  fair  results  for  the  long  wavalaagth  MStion  of  the 
fifth  mode.  8och  a  shall  theory  has  bean  prasaatad  by  liiraky'  “  for  a  siagla> 
layer  cylindrical  shall.  Resuming  tha  examination  of  Pigura  7  otts  aaos  that 
for  finita  wava  numbars  axial  motion  is  prasant  in  a  form  that  raaambles 
double >noda  thickoess  shear.  Howavar.  tha  axial  motiwi  shows  Uttls  tandancy 
toward  linearity,  indicating  that  (tha  aama  as  in  tha  kmar  modes)  a  shell  theory 
with  linear  axial  displacamant  variation  through  tha  thickaass  will  probably  not 
give  a  raasaaabla  approximation  for  larger  wava  numbers . 

In  summary  the  comparisons  of  tha  displacsmsaits  for  tha  four  lowest 
modes  shows  that  tha  sacoad  moda  yields  ths  poorest  comparison  between 
exact  aad  shall  theory.  Although  this  result  is  surprising,  It  could  have  bean 
anticipated  sines  tha  second  moda  disparaion  curve  davintos  more  from  that 
predicted  by  shall  theory  for  loarer  vaduas  of  A  than  for  tha  other  modes.  It 
would  appear  that  tha  incorporation  of  thickness -stretch  deformations  is  tha 
next  aasantial  step  in  improving  the  shall  theory  for  high  frequency,  short 
wavelength  use. 


-27- 


C.  STREMES 

FlfttT«a  8  through  11  daacrib*  tho  axiul  atroas  ab4  tha  ah«ar  atraaa 

di=trib::tic=: .  A*  dlapl^Mmanka  tha  ahaolttta  magnttuda  of  ^a 

atraaaaa  cannot  ba  datarminad.  Tha  atraaaaa  from  both  tha  aitact  and  ahaU 
thaonaa  wara  normaiiaad  ao  thai  um  girv^Uat  axial  stress  st  ths  enter  adgaa 
or  at  tha  Intarfaca  was  aat  aqual  to  aithar  al  .0»  This  also  fixas  tha  magaltttdaa 
of  tha  shaar  atraaaaa. 

Fifuraa  8  and  9  show,  raapactivaly.  tha  first  mods  axial  and  shaar 
strasaaa.  It  la  scan  in  Fig.  8  that  for  d  s  3  thara  is  good  agraatnant  with  tha 
shall  thaory.  Howavar,  as  d  ineraasas  tha  agraamant  datarioratas.  By  d  »  7 
thara  is  only  fair  agraamant,  and  by  d  =  9  thara  is  almost  nona.  Tha  lattar 
point  is  confirmed  by  tha  previous  obsarvatlcn  that  by  d  =  9  there  is  a  Rayleigh 
wave  formed  in  tha  outer  medium.  Figure  9  is  probably  of  more  interest  since 
an  axpactad  mode  of  failure  of  m\tltilayarad  sheila  is  in  shear  failure  of  tha 
bond,  Tha  shear  atraas  distrilution  follows  roughly  the  sama  pattam  as  the 
axial  stresses,  being  close  to  tha  shall  thaory  for  d  a  8  and  being  totally 
dilfaroBt  by  d  «  9. 

Figure  10  shows  tha  second  mode  strassa*.  As  might  ba  axpactad  thara 
is  not  as  good  an  agraamant  as  with  tha  first  mode.  By  d  <  3.0  tha  shall  thaory 
stress  distributions  are  beginning  to  diverge  aharply  front  tha  exact  ones. 

Figitva  11  ahowa  third  mode  atraaa  distributions.  Thera  is  good  agrss* 
mant  with  shall  thaory  only  for  d  s  1.  By  d  «  2  thara  ia  conaldarabla  divargance 
in  tha  atraaaaa,  and  by  d  «  3  tha  exact  theory  curves  bear  little  resemblance 
in  shapa  So  those  from  tha  shall  thaory.  Tha  fourth  mode  curves,  whlls  not 
prasentad,  show  s  pattam  similar  to  tha  third  mode.  It  is  worth  noting  that 
aithougb  ths  stress  comparison  plots  ara  ussful  for  illustrating  qualitative 
diifaroncaa  between  the  thaorica  tha  quantitative  tuts rp relation  can  ba  quite 
aanait*va  to  tha  particular  normalisation  convention  that  ia  adopted. 

In  concluaion  it  is  felt  that  if  any  shall  thaory  ia  to  ba  improved  by  tha  in¬ 
troduction  of  more  dependant  variables,  thicknaa a- stretch  motion  ia  aaaantial. 
For  atiU  higher  modes  to  yield  agraamant  it  would  ba  nacaaaary  to  add  additional 
•hear  deformations  corraapondlng  to  warping  of  initially  plana  cross  sections. 


•  SHEU  THEORY  - ELASTICITY  THEORY 

Fig.  8.  First  Mods  Axial  Strsaa  Oistribulions 


0  I  -I  0  I  -( 

SHEAR  STRESS 

SHELL  THEORY  - ELASTICITY  THEORY 


Fig.  9.  Firot  Mode  Stieat  Stress  Distributions 
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SHEAR  STRESS 


—  SHELL  THEORY  - EUSTICITY  THEORY 


Fig.  11.  Third  Mode  Streiis  Di»tribu lions 
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IV.  CONCLUSIONS 


The  rybsarvations  presented  h«r«  on  tb«  dynamics  of  o  two-layn*^ 
eyliudar  aro  spacifically  partinont  to  axially  symmotHc  motiOM  of  a  ralativaly 
tttia-'grallad  cylinder.  The  equetio^*.  K<M*N»v«r.  ara  In  no  way  raatrictad  to 
sholl-llko  goometry.  The  observations  are  also  pertinent  to  the  shell  theory 
of  Ref.  1  insofar  as  it  is  in  agreeniei^  with  the  exact  theory. 

In  many  respects  the  results  are  similar  to  those  for  a  homogeneous 
cylinder.  In  the  lowest  mode  very  long  waves  propagate  with  a  finite  phase 
velocity.  Infinite  wavelength  vibrations  in  the  second  mode  occur  with  a 
finite  cutoff  frequency  corresponding  to  exteneional  ring  vibrations  of  the 
cylinder.  Cutoff  frequencies  for  the  higher  modes  correspond  to  simple 
thickness* shear  or  thicknese-etretch  vibrations.  For  small  enough  wave¬ 
lengths  and/or  for  the  higher  modes  a  good  description  of  the  motion  is  obtained 
using  the  simpler  equations  for  the  propagation  of  straight  creeled  waves  in  a 
plans,  lediuin  of  the  same  thickness. 

However,  unlike  a  homogeneous  cylinder  a  two-layered  cylinder  does  not 
have  equal  surface  wave  velocitiee  at  its  inner  and  outer  eurfacae.  When  the 
flower  layer  is  much  thinner  this  leads  to  a  relative  maximum  at  intermediate 
wavelength  of  the  first  mode  phase  velocity  vs  wave  number  curve.  For  larger 
wave  numbers  the  phase  velocity  is  decreasing,  and  it  thus  approaches  from 
above  the  Rayleigh  wave  velocity  of  the  slower  medium.  This  charactsr  of  the 
phase  velocity  dispersion  relation  is  not  matched  even  qualitatively  by  a  Timo¬ 
shenko-type  shell  theory.  Presumably  duplication  of  this  bshavior  would  require 
a  higher  order  shell  theory  accommodating  cross  section  distortions  such  that 
initisklly  plane  cross  sections  of  either  layer  no  longer  remain  plane  in  the 
deformed  shell. 

Some  further  conclusions  pertinent  to  the  development  of  higher  order 
shell  theories  are  in  order.  The  displacement  and  etreee  diatributione  presented 
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here  reinforce  the  conclueion  that  good  roeolts  ar*  ohtainad  by  incorporatiag 

•  *  .4.  ^  s  t  a%«  >4  ^ 

aneer  aeiomi«iioo  UKUV^auauiy  in  v««;u  *sy«i  4-.  ^ 

it  appeare  that  if  oaa  deairod  to  oxtand  tha  raaga  of  applicability  e£  tha  Ihaei'f 
of  Ref.  1  by  adding  more  dependant  variablas  the  next  logical  etap  would  ba 
incorporation  of  thickoese- at  retch  daforonation  individually  in  each  l*yar. 
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C<|uaiif»iu  rssaltlaf  from  tlM  MtisfActloa  at  traa  MUtec* 
boundary  eoaditiona.  Eqa.  IS.  of^aikr  in  tbia  Appaadlx;  tkaif  atax 

Vii-® 

♦  ®2»‘^®20  tA-»> 

^l‘‘l‘’‘lO  ’  *  ^l^l^lO  '  10  •  ^**2‘^20  ■  •*•‘*'^*0 

tA-4) 

Aj<jl|0(«|a)  -  Bj<jKj<k«j*)  ♦  ‘  Dj*jKj^Sj») 

«  Aj^jljik^ja)  -  BjtjKjItfja)^  C.^ijIj^fcSja)  -  DjtjKjOtSja) 

(A-5) 

Ajlo(k*ja)  ♦  B,K,,0i«ja)  ♦  Cj6*I^(k6,a)  ♦  Dj6*K^(k*|a) 

.  Vo®'  V*  *  »2“o®'‘2**  *  *  D2®2Ko®‘*2*> 

^2**22  '*■  ®2°22  *  ^2*^22  ■*“  °2*22  *  ® 
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^^21  “l"**  '  ‘=2’»U  -  ®S*M  •  • 


IkmmnmMmm  •<« 


combtatttiOBB  of 


(A.9) 


"ij  •  (‘  *  •  **? 

•  **f  [‘o*‘l*'  ■  tI^] 
Sy  .»f[KofcV>*T^*] 

Ujj  •  2c^K|(k«jO} 

•  ‘i(‘  ♦  ‘f)‘i*V> 


i  toko*  OQ  valuoo  1  Obd  2.  Tho  oulMicrlpt  j  takoo  on  valuoo  1.  0.  and  2. 
wkleh  corroopooda  to  multiplying  tha  a^'gumanto  by  1  -f  X^,  1.  and  I  ~  X^, 
raopactivaly:  a.  g. . 
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Ha«a  Xj  •  (kj/a)  and  X^  ■  ^/a). 
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CYLINDER 
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Tho  Unaar  thaory  of  alaaticity  it  usad  to  invastif  ata  aaially  aymmatrlc  wava 
propagation  in  an  infiaitaly  long  taro>layarad  cylindar.  Eai^  matarial  ia 
takM  to  ba  homogaaaoua  and  iaotropic-  A  parfact  bond  ia  aaaomad  at  tha 
intarfaea,  arhila  dta  innar  and  outar  boundarioa  of  tha  compoaita  cylindar 
ara  traatad  aa  tr action -fraa.  Tha  diaparaion  datarmiaaat  ralatiag  phaaa 
valocity  and  wava  numbar  for  a  harmooic  train  of  aravaa  aatlafylnf  thaaa 
boundary  condittona  ia  praaantad.  Tha  charactar  of  tha  diaparaion  aquation 
i.  invaatigatad  analytically  and  numarically .  Straaa  and  diaplacamant 
diatributlona  ara  alao  praaantad  for  tha  numarical  anampla.  Compariaona 
ara  mada  with  an  appronimata  aolution  of  tha  aama  problam  obtainad  by 
maana  of  a  thin  ahall  thaory  incorporatiag  thiekaaaa-ahaar  daformatioit  of 
aach  layar. 
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